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Abstract 

The problem of electron-electron lifetime in a quantum dot is studied be- 
yond perturbation theory by mapping it onto the problem of localization in 
the Fock space. We identify two regimes, localized and delocalized, corre- 
sponding to quasiparticle spectral peaks of zero and finite width, respectively. 
In the localized regime, quasiparticle states are very close to single particle 
excitations. In the delocalized state, each eigenstate is a superposition of 
states with very different quasiparticle content. A transition between the two 
regimes occurs at the energy ~ A(^/ln g) 1 ^ 2 , where A is the one particle 
level spacing, and g is the dimensionless conductance. Near this energy there 
is a broad critical region in which the states are multifractal, and are not 

described by the Golden Rule. 
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Quasiparticle in a Fermi-liquid is not an eigenstate: it decays into two quasiparticles 
and a hole. In an infinite clean system, by using the Golden Rule (GR), quasiparticle decay 
rate is estimated as 7(e) ~ e 2 /ep, where e is quasiparticle energy and ep is Fermi energy 
|Ij]. However, in a finite system all single- and many-particle states are discrete. In this 
case, quasiparticles may be viewed as wave-packets constructed of such states, the packet 
width being determined by the quasiparticle lifetime in an infinite system: 8e ~ 7(e). In this 
paper we attempt at making the relation between quasiparticles and many-particle states 
more precise, and find that at different energies it has different meaning. 

Conventionally, quasiparticles are well defined provided 7(e) <C e. However, to resolve 
quasiparticles in a mesoscopic system, a more stringent condition is required: 7(e) < A, 
the quasiparticle level spacing. As an example, we consider quasiparticle peaks in tunneling 
conductance of a quantum dot 0|§. The peaks are observed in non-linear conductance at 
certain bias, and are interpreted in terms of energy dependence of the quasiparticle tunneling 
density of states (DOS), so that each peak corresponds to a "quasiparticle state," and its 
width measures the lifetime of the state. Below, in our discussion of the width in terms of 
constituting eigenstates, we ignore any contributions to the quasiparticle decay due to finite 
escape rate, phonons, etc. and consider an isolated Fermi-liquid. 

The meaning of quasiparticle lifetime needs clarification: strictly speaking, since quan- 
tum dot is a finite system, any many-particle eigenstate gives rise to an infinitely narrow 
conductance peak. However, we will see that only a small fraction of those states have sig- 
nificant overlap with one-particle excitations, and thus can be detected by a finite sensitivity 
measurement. Under certain conditions, these strong peaks group into clusters of the width 
~ 7(e), that can be interpreted as quasiparticle peaks. 

Before discussing possible regimes let us review the GR approach. Recently Sivan et al. 
0, adopting the quasi-particle picture to a finite size geometry and relying on the earlier 
work || on electron-electron scattering rate in diffusive conductors found that 

7(e) » A • (e/gA) 2 , e < gA , (1) 
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where A is the mean single-particle level spacing near Fermi level and g 3> 1 is the dimen- 
sionless conductance, for a finite system defined by g = E c /A, where E c is the Thouless 
energy (inverse time of diffusion through the system). The decay rate ([I]) is much larger 
than in a clean Fermi-liquid, however, at e < E c one has 7(e) <C A, implying that the 
quasiparticle states can be resolved. 

However, the GR can be used to evaluate lifetime only when the density of final states 
is sufficiently large, so that the GR decay rate is larger than the level spacing of final 
unperturbed states. Otherwise, the GR will not give the decay rate, but rather just a first- 
order perturbation correction to the energy of a given eigenstate. It is important to realize 
that in our problem, since a quasiparticle decays into three quasiparticles, the density of 
relevant final states, z/ 3 (e) = e 2 /(2A 3 ), is much smaller than that of all many-body states. 
The interaction matrix element V in the GR leading to Eq. [1] is of the order of A/g (see 
below), which should be compared to the three-particle level spacing l/z/3. Therefore, the 
GR is not applicable unless e > e* = A^fg. Note that, since e* A, there are many states 
whose lifetime is not given by GR. 

At e e*, when matrix elements are much smaller than the spacing l/z/3, the quasipar- 
ticle states do not decay: they are just slightly perturbed one-particle states. Hence they 
produce strong and very narrow conductance peaks, some of which may have weak satellites 
due to coupling to many-particle states. As e approaches e* from below, the number of 
the satellites rapidly increases, and at e ^> e*, they altogether form finite width peaks, well 
described by the GR. 

For a quantitative description of the whole interval < e < E c (including the vicinity 
of e*), it is both interesting and instructive to explore the analogy of this problem with the 
Anderson localization. This will be the goal of this paper. 

Extension of the traditional localization problem to few interacting particles has received 
much attention recently. The study of the two-particle case, started by Dorokhov J7|, was 
further advanced by Shepelyansky Imry || and Pichard et al. [|10|| , with extensions to 
more particles. Some of the energy scales encountered below, such as V and l/z/3, have 



already been discussed in the context of those papers. In contrast, here we deal with the 
states extended throughout the whole finite system. The localization we consider occurs in 
the Fock space of many-body states, rather than in the real space. 

Distance in Fock space is a concept we introduce to measure closeness of many-body states. 
Consider a generic two-body interaction Hamiltonian in a secondary quantized form 

H + Hi = e *c+c a + V^f c ^ c 5 c /3C a - (2) 

a a/375 

The many-body problem is formulated in the Fock space, by choosing as a basis {^n} - 
the Slater determinants constructed out of the iV-particle Fermi vacuum \N): 

*N = c+ 2m ...c+ m+l c am ...c ai \N) . (3) 

Any state \I/ n can be represented as a string with entries 1 and labeling the single particle 
states which do or do not participate in v&n, respectively. Let us define the distance between 
two states and ^/' N as the number of positions in which the corresponding strings differ. 

Thus defined, the distance in the Fock space can be used to map the lifetime problem 
to an appropriate Anderson localization problem. It is useful to think of the states \l/jv, 
the eigenstates of H , as of "site orbitals" in the Fock space, each having an on-site energy 
e <*2m + ••• + e a m+1 + \ca m \ + ••• + l e aj- These sites are interconnected by the interaction Hi, 
which we think of as hopping in the Anderson problem (the diagonal part of Hi is added to 
Ho by using the Hartree-Fock method). The point is that the two-body interaction matrix 
element (\E'Ar|7ii| , I'^ v ) is non-zero only if the distance between the states and ^/' N equals 
0, 2 or 4. We construct a network in the Fock space by connecting all orbitals ^jy which 
are at a distance 2 from each other. 

Below we study the localization problem on this network. We show that the energies 
at which the GR is relevant correspond to states extended over the network, whereas at 
lower energies the hopping does not form enough resonances between neighboring orbitals, 
resulting in localization of eigenstates near original sites. Unlike in other localization prob- 
lems, as the energy is lowered, localization develops very gradually, within a large critical 
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region, e** < e < e*, where e** ~ e*/y/\ng and A < e** < e* < £7 C . This is a result of 
the particular structure of the network in the Fock space. Above the localization threshold 
(and beyond the critical region) we recover the GR picture |J with finite width quasiparticle 
conductance peaks. Each of them is formed by a huge number of eigenstate peaks, and even 
a small external broadening will smear this fine structure. In the localized regime e < e** 
the single-particle density of states will consist of isolated 5— function peaks. In the critical 
region, e** < e < e*, the peaks have irregular structure which turns into simple Lorentzian 
above e*. 

The meaning of localization in the Fock space is that a localized state is practically 
identical to a single-particle excitation (or a superposition of very few quasiparticle states). 
The energy of each of the constituent quasiparticles represents a good quantum number, 
whereas for the extended delocalized states only total energy is conserved. The transition 
is of the Anderson type because the two-body Hamiltonian is local in the Fock space: it 
couples only the orbitals of similar quasiparticle content. 

It is worth remarking that the hierarchical organization of many-particle states proved 
to be a useful picture in studies of compound nuclear reaction rates flH| |. The "doorway 



states" introduced in the nuclear reaction studies, although serve a different purpose, are 
related to our network construction. 

Hopping over the network in the Fock space: The two-body interaction matrix elements in 
(§) are given by 

V;f = JJdxdx'V(x - x'WsixW^M^Mx') • (4) 

To evaluate the matrix elements, let us consider diffusive disorder and a short range inter- 
action, V(x — x') = XAV5(x — x'), where V is the volume, and A ~ 1 is the dimensionless 
interaction strength. For a ^ (3 ^ 7 7^ 5, V^f is a random quantity with zero average. The 



root-mean-square V may be evaluated [12], e.g., by using the diagram shown in Fig. l.a. 



In the absence of time reversal symmetry we obtain 
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A 2 2 '-y 2 

^ = « = ^E^. (5) 

where 7 m are eigenvalues of the diffusion operator. (By definition, E c — 71.) In deriving 
(|) we assume small single-particle energies: e a (i3,-y,s) < ^c- The magnitude of V decreases 
algebraically when the differences between the single-particle energies exceed E c ; below we 
shall ignore such contributions. 

The network in the Fock space organizes all states in a hierarchy. Let \N — 1) be the 
ground state of N — f particles. The states T a = c+ \N — 1) representing one particle 
added in the state a form the first generation of the hierarchy (see Fig. l.b). The states 
T"^ = c^CpC^N — 1), representing two particles and one hole, form the generation 3. 
Similarly, the generation 5 is formed by T"^ 7 = c^CpC^C\c^\N — 1), etc. The two-body 
interaction 7i\ couples only the states of near generations, so that any given state from 
generation 2n — 1 is connected only to states from generations 2n + 1, 2n — 1, or 2n — 3. 
This implies that connected states are a distance 2 from each other. 

Consider now a state of generation 1, with an on-site energy e. The DOS in generation 
3 accessible by "hopping," having the same energy is ^(e) = e 2 /(2A 3 ). For higher-order 
generations (more excited quasi-particles) , the DOS rapidly increases: for the (2n + l) st 



generation (where n < n max ps ye/A) it goes as (e/A) 2n /(2n)!. However, we should focus 
only on those states of generation (2n + 1) which are directly accessible from a given state of 
generation (2n — 1). The density of such states is much smaller, and is given by z^n.+i = 
We note that from a state in the generation (2n — 1) it is also possible to hop to some states 
of the same generation, and to some states of the previous generation (2n — 3). Respectively, 



DOS associated with these processes is yz/ 3 /A and n(n — l)(2n — 3)/A. For n n max the 
number of such hopping processes is parametrically smaller than the number of states in 
the next generation accessible by hopping. We thus obtain a picture which is quite close to 
that of a Cayley tree: each "site" of the (2n — l) st generation branches out to K n sites of 
the next generation. (The number of couplings to the sites of the same or of the previous 



generations is much smaller, and thus can be ignored fl3| .) The branching number is given 



by integrating the effective DOS over the energy interval E c where the hopping parameter 
V is energy independent. We obtain the branching number 

K n « g 3 /(Qn) . (6) 

Note that K n depends on the hierarchy level of the tree. 

Decreasing the branching number K n with increasing n makes the network effectively 
finite. In order to simplify the consideration, below we will consider an infinite Cayley 
tree with constant branching number K — K\ — g 3 /Q. More realistic model, taking into 
account finite size of the tree and n-dependence of the branching number, will be considered 



elsewhere ||14j| . Qualitatively, as a result of the finiteness of the network, the localization 
transition we find in the infinite network will be smeared to a crossover. Up to the finite 
size smearing effects, all the results obtained for the infinite Cayley tree will remain intact. 

The model we are interested in was solved by Abou-Chacra, Anderson and Thouless 
I5| . They considered localization on a Cayley tree with the on-site energies from a uniform 



distribution in the interval [— W, W], and constant hopping amplitude V. By studying 
fix-points of the mapping of self-energies computed recursively using the hierarchy of the 
Cayley tree, it was found that delocalization occurs at V ~ W/(KlnK), where K is the 
tree branching number. This result is in apparent disagreement with the Anderson criterion 
of localization, V ~ W/ K, obtained by comparing the hopping amplitude with the spacing 
of the on-site energies in the nearest neighbor shell. 



To understand the significance of linK, instead of following the exact solution ||15|| , let 
us study the statistics of resonances of the tree sites appearing due to hopping. Starting at 
a site of the first generation, the amplitude (at energy e) connecting this site to a given site 
at the (2n + l) s * generation, in lowest order in V is given by 

n y 

1=1 fc fc * 

Let us consider the probability, p(n,C), that this amplitude is significant, namely that 
|yi n | exceeds given finite C. We define In \A n \ = nln^ + Y n , where Y n = Y% =1 yi and 



Vi = l n 1 I • Assuming that the relevant on-site energies are uniformly distributed in the 
interval — W < e — e« < W, the probability distribution of yi is P{yi) = exp{— where 
< n« < oo. Fourier transforming P, taking the n t/l power and Fourier transforming back, 
we obtain the distribution function of Y n , (0 < Y n < oo). Since \A n \ = (-^)™ exp{F„}, the 



distribution function of |A„| will be 

(V/W) n 



P(\An\) 



in -IV. IAJ 2 



n-1 



(8) 



Thus the probability p(n, C) = dA P(A) for W ^> V, C ~ 1 is given by 



p{n,C) 



1 
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(9) 



(n-i)! cm (c(f; 

The probability that none of the trajectories connecting a site in the first generation to 
sites in the generation (2n + 1) carries a large amplitude is given by 



(l-p(n,C)) K =exp(-/ B ), 

where for p(n, C) < 1, /„ « K n p(n, C). From Eq. ^ for n ^> 1 one obtains 

1 1 



(10) 
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If f n increases at large n, then, eventually, at higher generations one has f n 3> 1, i.e., strong 
coupling to generation 1. A transition to the localized phase takes place when the expression 
in the square brackets in Eq. |TT| assumes the value 1, which gives the criterion quoted above. 

To apply this result to our problem, we replace V"f by V and approximate the density 
of states in the generation (2n + 1), accessible from a state (of energy e) in the generation 
(2n — 1), to be uniform in the interval [e - W, e + W], V < W < e, and equal to Kj (2W) = 



z/ 3 (e). We then find that the transition occurs at the energy (16 

e** ss (\b d 2e)- 1 ' 2 



In g 



(12) 



At energies just above e** the first generation is not very well connected with the next 
few generations (f n < 1 for small values of n). The condition for all generations to be 
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well connected with the generation 1 is fi = KV/{WC) > 1. This coincides with a naive 
implementation of the Thouless criterion ||17|| , yielding the second energy scale 



Let us discuss the meaning of the various regimes. In the localized phase (e < e**) the 
first generation is weakly connected with the rest of the network. Therefore, at such energies 
the exact many-body states are close to Slater determinants, or to superpositions of few 
determinants. A mathematical description of a single particle injection into a dot involves 
projecting a single-particle state onto exact eigenstates of the system. In the localized phase 
each single-particle state will have a significant overlap with one (or few) exact eigenstates, 
producing a few resolved 5-function peaks in the spectrum of the single-particle DOS. At 
e > e* all generations are well connected. Due to the huge density of multiparticle states, the 
states of generation 1 can be thought of as being effectively well coupled to the continuum. 
This justifies the GR result, Eq. [Tl in this energy range [0,0. Each single particle peak 
associated with generation 1 is replaced by a cluster of a large number of many-particle peaks, 
altogether forming a Lorentzian envelope. For e < E c , the width of the envelope is less than 
A, and thus the "quasiparticle states" can be resolved in, e.g., transport measurements ||. 
For intermediate energies e** < e < e* there are still many peaks in a cluster. However the 
probability that a particular generation is represented in a given cluster is small. As a result 
the widths of clusters as well as the shapes of their envelopes will strongly fluctuate from 
peak to peak. 

Following from our discussion, there are some interesting implications to the relation of 
Quantum Chaos and Anderson Localization. As long as quasiparticle peaks can be resolved, 
the many-particle spectrum is not trully chaotic in the Wigner-Dyson sense. (For example, 
from our analysis, different quasiparticle states "do not talk to each other.") Another 




(13) 



A convenient expression for f n is, then, 




(14) 
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observation is that an extended state in the many-body problem is not necessarily ergodic, 
but can be extended over a small fraction of the whole space (in our case, over a subtree of 
the Cayley tree). 
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FIG. 1. (a) The diagrams used to evaluate the mean square of the matrix element in Eq. || . 
(b) Schematic representation of the Cayley tree in the Fock space of many-body states. Different 
generations are shown. 
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